Analysis of changes occurring in Codon Positions due to mutations through the cellular automata transition rules (preprint)

Variation in the nucleotides of a codon may cause variations in the evolutionary patterns of a DNA or amino acid sequence. To address the capability of each position of a codon to have non-synonymous mutations, the concept of degree of mutation has been introduced. The degree of mutation of a particular position of codon defines the number of non-synonymous mutations occurring for the substitution of nucleotides at each position of a codon, when other two positions of that codon remain unaltered. A Cellular Automaton (CA), is used as a tool to model the mutations of any one of the four DNA bases A, C, T and G at a time where the DNA bases correspond to the states of the CA cells. Point mutation (substitution type) of a codon which characterizes changes in the amino acids, have been associated with local transition rules of a CA. Though there can be [Formula] transitions of a 4-state CA with 3-neighbourhood cells, here it has been possible to represent all possible point mutations of a codon in terms of combinations of 16 local transition functions of the CA. Further these rules are divided into 4 classes of equivalence. Also, according to the nature of mutations, the 16 local CA rules of substitutions are classified into 3 sets namely, No Mutation, Transition and Transversion. The experiment has been carried out with three sets of single nucleotide variations(SNVs) of three different viruses but the symptoms of the diseases caused by them are to some extent similar to each other. They are SARS-CoV-1, SARS-CoV-2 and H1N1 Type A viruses. The aim is to understand the impact of nucleotide substitutions at different positions of a codon with respect to a particular disease phenotype.

Non-Uniform Aspects of SARS-CoV-2 Intraspecies Evolution Reopen Questions on Its Origin (preprint)

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no hypothesis has managed to identify the origin, and the issue has resurfaced. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins across different continents comprising 24 geo-locations. The results showed an evenly uneven distribution of unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across the 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations we have considered. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and to be prepared to meet the challenges of potential future pandemics.

An Issue of Concern: Unique Truncated ORF8 Protein Variants of SARS-CoV-2 (preprint)

Open reading frame 8 (ORF8) protein is one of the most evolving accessory proteins in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits presentation of viral antigens by the major histocompatibility complex class I (MHC-I) and interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 to evade immunity and replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein defines the B.1.1.7 lineage of SARS-CoV-2, which is engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) due to the Q27STOP mutations were identified among 49055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents which includes Africa, Asia, Europe and South America. Based on various quantitative features such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, a collection of nine possible T-ORF8 unique variants were defined. The question of whether T-ORF8 variants work similarly to ORF8 has yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Implications Derived from S-Protein Variants of SARS-CoV-2 from Six Continents (preprint)

Spike (S) proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are critical determinants of the infectivity and antigenicity of the virus. Several mutations in the spike protein of SARS-CoV-2 have already been detected, and their effect in immune system evasion and enhanced transmission as a cause of increased morbidity and mortality are being investigated. From pathogenic and epidemiological perspectives, spike proteins are of prime interest to researchers. This study focused on the unique variants of S proteins from six continents Asia, Africa, Europe, Oceania, South America, and North America. In comparison to the other five continents, Africa (29.065%) had the highest percentage of unique S proteins. Notably, only North America had 87% (14046) of the total (16143) specific S proteins available in the NCBI database(across all continents). Based on the amino acid frequency distributions in the S protein variants from all the continents, the phylogenetic relationship implies that unique S proteins from North America were significantly different from those of the other five continents. Overtime, the unique variants originating from North America are most likely to spread to the other geographic locations through international travel or naturally by emerging mutations. Hence it is suggested that restriction of international travel should be considered, and massive vaccination as an utmost measure to combat the spread of COVID-19 pandemic. It is also further suggested that the efficacy of existing vaccines and future vaccine development must be reviewed with careful scrutiny, and if needed, further re-engineered based on requirements dictated by new emerging S protein variants.

Carbon-Based Nanomaterials: Promising Antiviral Agents to Combat COVID-19 in the Microbial Resistant Era (preprint)

Therapeutic options for the highly pathogenic human Severe Acute Respiratory Syndrome-Coronavirus 2 (SARS-CoV-2) causing the current pandemic Coronavirus disease (COVID-19) are urgently needed. COVID-19 is associated with viral pneumonia and acute respiratory distress syndrome causing significant morbidity and mortality. The proposed treatments for COVID-19, such as hydroxychloroquine, remdesivir and lopinavir/ritonavir, have shown little or no effect in the clinic. Additionally, bacterial and fungal pathogens contribute to the SARS-CoV-2 mediated pneumonia disease complex. The antibiotic resistance in pneumonia treatment is increasing at an alarming rate. Therefore, carbon-based nanomaterials (CBNs), such as fullerene, carbon dots, graphene, and their derivatives constitute a promising alternative due to their wide-spectrum antimicrobial activity, biocompatibility, biodegradability and capacity to induce tissue regeneration. Furthermore, the antimicrobial mode of action is mainly physical (e.g. membrane distortion), which is characterized by a low risk of antimicrobial resistance. In this review, we evaluated the literature on the antiviral activity and broad-spectrum antimicrobial properties of CBNs. CBNs had antiviral activity against 12 enveloped positive-sense single-stranded RNA viruses similar to SARS-CoV-2. CBNs with low or no toxicity to the humans are promising therapeutics against COVID-19 pneumonia complex with other viruses, bacteria and fungi, including those that are multidrug-resistant.

Clade GR and Clade GH Isolates in Asia Show Highest Amount of SNPs (preprint)

Clades are monophyletic groups composed of a common ancestor and all its lineal descendants. As the propensity of virulence of a disease depends upon the type of clade the virus belongs to and it causes different fatality rates of disease in different countries, so the clade-wise analysis of SARS-CoV-2 isolates collected from different countries can illuminate the actual evolutionary relationships between them. In this study, 1566 SARS-CoV-2 genome sequences across ten Asian countries are collected, clustered, and characterized based on the clade they belong to. The isolates are compared to the Wuhan reference sequence (Accession no:MN996528.1) to identify the mutations that occurred at different protein regions. Structural changes in amino acids due to mutations lead to functional instability of the proteins. Detailed clade-wise functional assessments are carried out to quantify the stability and vulnerability of the mutations occurring in SARS-CoV-2 genomes which can shade light on personalized prevention and treatment of the disease and encourage towards the invention of clade-specific vaccines.

Variability of Accessory Proteins Rules the SARS-CoV-2 Pathogenicity (preprint)

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) which is pandemic with an estimated fatality rate less than 1% is ongoing. SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 with putative functions to manipulate host immune mechanisms such as interferons, immune signaling receptor NLRP3 (NOD-, LRR-, and pyrin domain-containing 3) inflammasome, inflammatory cytokines such as interleukin {beta} (IL-1{beta}) are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins of all complete proteomes (available as of October 26, 2020, in the National Center for Biotechnology Information depository) of SARS-CoV-2, were observed across six continents. Across all continents, the decreasing order of percentage of unique variations in the accessory proteins was found to be ORF3a>ORF8>ORF7a>ORF6>ORF10>ORF7b. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. This finding suggests that the wide variations of accessory proteins seem to govern the pathogenicity of SARS-CoV-2, and consequently, certain propositions and recommendations can be made in the public interest.

Forecasting COVID-19 cases at the Amazon region: a comparison of classical and machine learning models (preprint)

BACKGROUNDSince the first reports of COVID-19, decision-makers have been using traditional epidemiological models to predict the days to come. However, the enhancement of computational power, the demand for adaptable predictive frameworks, the short past of the disease, and uncertainties related to input data and prediction rules, also make other classical and machine learning techniques viable options. OBJECTIVEThis study investigates the efficiency of six models in forecasting COVID-19 confirmed cases with 17 days ahead. We compare the models autoregressive integrated moving average (ARIMA), Holt-Winters, support vector regression (SVR), k-nearest neighbors regressor (KNN), random trees regressor (RTR), seasonal linear regression with change-points (Prophet), and simple logistic regression (SLR). MATERIAL AND METHODSWe implement the models to data provided by the health surveillance secretary of Amapaa, a Brazilian state fully carved in the Amazon rainforest, which has been experiencing high infection rates. We evaluate the models according to their capacity to forecast in different historical scenarios of the COVID-19 progression, such as exponential increases, sudden decreases, and stability periods of daily cases. To do so, we use a rolling forward splitting approach for out-of-sample validation. We employ the metrics RMSE, R-squared, and sMAPE in evaluating the model in different cross-validation sections. FINDINGSAll models outperform SLG, especially Holt-Winters, that performs satisfactorily in all scenarios. SVR and ARIMA have better performances in isolated scenarios. To implement the comparisons, we have created a web application, which is available online. CONCLUSIONThis work represents an effort to assist the decision-makers of Amapa in future decisions to come, especially under scenarios of sudden variations in the number of confirmed cases of Amapa, which would be caused, for instance, by new contamination waves or vaccination. It is also an attempt to highlight alternative models that could be used in future epidemics.

Possible transmission flow of SARS-CoV-2 based on ACE2 features (preprint)

Angiotensin-converting enzyme 2 (ACE2) is the cellular receptor for the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that is engendering the severe coronavirus disease 2019 (COVID-19) pandemic. The spike (S) protein receptor-binding domain (RBD) of SARS-CoV-2 binds to the three sub-domains viz. amino acids (aa) 22-42, aa 79-84, and aa 330-393 of ACE2 on human cells to initiate entry. It was reported earlier that the receptor utilization capacity of ACE2 proteins from different species, such as cats, chimpanzees, dogs, and cattle, are different. A comprehensive analysis of ACE2 receptors of nineteen species was carried out in this study, and the findings propose a possible SARS-CoV-2 transmission flow across these nineteen species.

Notable sequence homology of the ORF10 protein introspects the architecture of SARS-COV-2 (preprint)

The global public health is endangered due to COVID-19 pandemic, which is caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2). Despite having similar pathology to MERS and SARS-CoV, the infection fatality rate of SARS-CoV-2 is likely lower than 1%. SARS-CoV-2 has been reported to be uniquely characterized by the accessory protein ORF10, which contains eleven cytotoxic T lymphocyte (CTL) epitopes of nine amino acids length each, across various human leukocyte antigen (HLA) subtypes. In this study, all missense mutations found in sequence databases were examined across twnety-two unique SARS-CoV-2 ORF10 variants that could possibly alter viral pathogenicity. Some of these mutations decrease the stability of ORF10, e.g. I4L and V6I were found in the MoRF region of ORF10 which may also possibly contribute to Intrinsic protein disorder. Furthermore, a physicochemical and structural comparative analysis was carried out on SARS-CoV-2 and Pangolin-CoV ORF10 proteins, which share 97.37% amino acid homology. The high degree of physicochemical and structural similarity of ORF10 proteins of SARS-CoV-2 and Pangolin-CoV open questions about the architecture of SARS-CoV-2 due to the disagreement of these two ORF10 proteins over their sub-structure (loop/coil region), solubility, antigenicity and change from the strand to coil at amino acid position 26, where tyrosine is present. Altogether, SARS-CoV-2 ORF10 is a promising pharmaceutical target and a protein which should be monitored for changes which correlate to change pathogenesis and clinical course of COVID-19 infection.

Pathogenetic Perspective of Missense Mutations of ORF3a Protein of SARS-CoV2 (preprint)

One of the most important proteins for COVID-19 pathogenesis in SARS-CoV2 is the ORF3a protein which is the largest accessory protein among others accessory proteins coded by coronavirus genome. The major roles of the protein include virulence, infectivity, ion channel activity, morphogenesis and virus release. The coronavirus, SARS-CoV2 is continuously evolving naturally and thereby the encoded proteins are also mutating rapidly. Therefore, critical study of mutations in ORF3a is certainty important from the pathogenetic perspective. Here, a sum of 175 various non-synonymous mutations in the ORF3a protein of SARS-CoV2 are identified and their corresponding effects in structural stability and functions of the protein ORF3a are studied. Broadly three different classes of mutations, such as neutral, disease and mixed (neutral and disease) type mutations were observed. Consecutive mutations in some ORF3a proteins are established based on timeline of detection of mutations. Considering the amino acid compositions over the ORF3a primary protein sequences, twenty clusters are detected based on K-means clustering method. Our findings on 175 novel mutations of ORF3a proteins will extend our knowledge of ORF3a, a vital accessory protein in SARS-CoV2, which would assist to enlighten on the pathogenicity of this life-threatening COVID-19.

SARS-CoV2 Envelope Protein: Non-Synonymous Mutations and Its Consequences (preprint)

In the NCBI database, as on June 6, 2020, total number of available complete genome sequences of SARS-CoV2 across the world is 3617. The envelope protein of SARS-CoV2 possesses several non-synonymous mutations over the transmembrane domain and (C)-terminus in 0.414\% of these 3617 genomes. The C-terminus motif DLLV has been changed to DFLV and YLLV in the proteins QJR88103 (Australia: Victoria) and QKI36831 (China: Guangzhou) respectively, which might affect the binding of this motif with the host protein PALS1.

Molecular conservation and Differential mutation on ORF3a gene in Indian SARS-CoV2 genomes (preprint)

A global emergency due to the COVID-19 pandemic demands various studies related to genes and genomes of the SARS-CoV2. Among other important proteins, the role of accessory proteins are of immense importance in replication, regulation of infections of the coronavirus in the hosts. The largest accessory proteins in the SARS-CoV2 genome is ORF3a which modulates the host response to the virus infection and consequently it plays an important role in pathogenesis. In this study, an attempt is made to decipher the conservation of nucleotides, dimers, codons and amino acids in the ORF3a genes across thirty two genomes of Indian patients. ORF3a gene possesses single and double point mutations in Indian SARS-CoV2 genomes suggesting the change of SARS-CoV2s virulence property in Indian patients. We find that the parental origin of the ORF3a gene over the genomes of SARS-CoV2 and Pangolin-CoV is same from the phylogenetic analysis based on conservations of nucleotides and so on. This study highlights the accumulation of mutation on ORF3a in Indian SARS-CoV2 genomes which may provide the designing therapeutic approach against SARS-CoV2.

Some Findings on Genes over SARS-CoV2 Genomes (preprint)

Coronaviruses are a large family of RNA viruses which cause respiratory infections ranging from the common cold to more severe diseases such as Middle East Respiratory Syndrome (MERS), Severe Acute Respiratory Syndrome (SARS) and COVID-19. This article highlights some key findings based on a thorough scanning of genes of 475 SARS-CoV2 genomes, including the co-presence of ORF7a and ORF8 over the 256 SARS-CoV2 genomes and the absence of the gene ORF7b over the 219 SARS-CoV2 genomes collected from various countries including India. The presence of the gene ORF7b is found in the SARS-CoV2 genomes containing the L-type strain which is reported to having much higher virulence as compared to the S-type strain.

On spatial molecular arrangements of SARS-CoV2 genomes of Indian patients (preprint)

A pandemic caused by the SARS-CoV2 is being experienced by the whole world since December, 2019. A thorough understanding beyond just sequential similarities among the protein coding genes of SARS-CoV2 is important in order to differentiate or relate to the other known CoVs of the same genus. In this study, we compare three genomes namely MT012098 (India-Kerala), MT050493 (India-Kerala), MT358637 (India-Gujrat) from India with NC_045512 (China-Wuhan) to view the spatial as well as molecular arrangements of nucleotide bases of all the genes embedded in these four genomes. Based on different features extracted for each gene embedded in these genomes, corresponding phylogenetic relationships have been built up. Differences in phylogenetic tree arrangement with individual gene suggest that three genomes of Indian origin have come from three different origins or the evolution of viral genome is very fast process. This study would also help to understand the virulence factors, disease pathogenicity, origin and transmission of the SARS-CoV2.